Fluctuation-Driven Molecular Transport in an Asymmetric Membrane Channel 
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Channel proteins, that selectively conduct molecules across cell membranes, often exhibit an 
asymmetric structure. By means of a stochastic model, we argue that channel asymmetry in the 
presence of non-equilibrium fluctuations, fueled by the cell's metabolism as observed recently, can 
dramatically influence the transport through such channels by a ratchet-like mechanism. For an 
aquaglyceroporin that conducts water and glycerol we show that a previously determined asymmetric 
glycerol potential leads to enhanced inward transport of glycerol, but for unfavorably high glycerol 
concentrations also to enhanced outward transport that protects a cell against poisoning. 

PACS numbers: 87.16.Uv, 05.40.-a, 05.10.Gg 



Living cells interact with their extracellular environ- 
ment through the cell membrane, which acts as a pro- 
tective permeability barrier for preserving the internal 
integrity of the cell. However, cell metabolism requires 
controlled molecular transport across the cell membrane, 
a function that is fulfilled by a wide variety of transmem- 
brane proteins, acting as passive and active transporters 
^J. Channel proteins as passive transporters facilitate 
the diffusion of specific molecules across the membrane 
down a free energy gradient. Active transporters con- 
duct molecules along or against the free energy gradient 
consuming for that purpose external energy. However, 
the plasma membrane of living cells is subject to a vari- 
ety of non-equilibrium, i.e., non-thermal processes, e.g., 
interaction with the cytoskeleton and with active mem- 
brane proteins 0,0. In this Letter we want to argue 
that in an active plasma membrane even passive chan- 
nel proteins can act as active transporters by consuming 
energy from non-equilibrium fluctuations fueled by cell 
metabolism 

The Escherichia coli glycerol uptake facilitator (GlpF) 
is an aquaglycerol channel protein, which transports both 
water and glycerol molecules, but excluding charged so- 
lutes, e.g., protons The recently determined 3D struc- 
ture of GlpF at atomic resolution 6] (Fig. [IJ has pro- 
vided much insight into the underlying microscopic mech- 
anism of molecular transport and selectivity through 
]. In particular, molecular dynamics (MD) stud- 
, [HI E| established that water and glycerol dif- 
fusion through GlpF is single file, and for biologically 
relevant periplasmic glycerol concentrations correlation 
effects between consecutive glycerol molecules are neg- 
ligible due to their large spatial separation. The corre- 
sponding potential of mean force (PMF) that guides 
the transport of glycerol through the channel is highly 
asymmetric reflecting the atomic structure of GlpF (see 
Fig. HJ, with a prominent potential well at the exter- 
nal (periplasmic) side and a constriction region with sev- 




eral pronounced potential barriers towards the internal 
(cytoplasmic) side of the channel jy, [l2[. Besides GlpF, 



cytoplasm 




0.2 
periplasm 



0.4 ,0.6 
x/L 



FIG. 1: (color online). Section through the glycerol con- 
duction pathway in GlpF and the corresponding asymmetric 
PMF (solid curve) reported in Q. 



there are several other porins which exhibit similar spa- 
tial asymmetry [ill LL3 > an d in spite of recent efforts, 
no biological function could be attributed to the asym- 
metry . Here we propose and demonstrate that under 
realistic physiological conditions, the asymmetry of GlpF 
furnishes active glycerol transport through a ratchet-like 
mechanism. In general, the ratchet effect refers to the 
generation of directed motion of Brownian particles in 
a spatially periodic and asymmetric (ratchet) potential 
in the presence of non- equilibrium fluctuations and/or 
externally applied time-periodic force with zero mean 
[lTl H^ . Since cell membranes are subject to non- 
equilibrium fluctuations, which span a wide range of time 
and length scales 0, , one expects a ratchet effect con- 
tribution to the transport through asymmetric channel 
proteins, such as GlpF, even if the PMF is non-periodic. 
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Assuming simplified, heuristic models, the ratchet effect 
has been invoked before to explain the functioning of ac- 
tive biomolecules, e.g., motor proteins and ATP hydrol- 
ysis driven ion pumps 0, |2fl|. To our knowledge this 
Letter uses for the first time a realistic, microscopically 
determined PMF to investigate the precise role of non- 
equilibrium fluctuations in facilitated transport. We find 
that as a result of channel asymmetry glycerol uptake, 
driven by a concentration gradient, is enhanced signif- 
icantly in the presence of non-equilibrium fluctuations. 
Furthermore, the ratchet effect-caused enhancement is 
larger for the outward, i.e., from the cytoplasm to the 
periplasm, flux than for the inward one, suggesting that 
non-equilibrium fluctuations also play an important role 
in protecting the interior of the cell against excess uptake 
of glycerol. 

Glycerol transport through GlpF can be modeled in 
terms of overdamped Brownian motion along the axis 
of the channel as a result of the concentration gradi- 
ent established at the ends of the channel. The inter- 
action of a diffusing glycerol molecule with the protein, 
solvent, lipid and other glycerol molecules is taken into 
account through the PMF, U(x), as determined from 
steered molecular dynamics simulations by employing the 
Jarzynski equality |S| . The motion of a glycerol molecule 
inside GlpF and in the presence of an external force F(t) 
in the strong friction limit is described by the Langevin 
equation (LE): jx = f(x) + + F(t), where 7 is the 
friction coefficient, f(x) = —U'{x) is the deterministic 
force derived from the PMF, and £(£) is the Langevin 
force due to the equilibrium thermal fluctuations. As 
usual, £(t) is a Gaussian white noise with (£(£)) = and 
(f (0)) = 2D^ 2 5(t), where 5(b) is the Dirac-delta func- 
tion, and D is the effective diffusion coefficient of a glyc- 
erol molecule inside GlpF. According to the fluctuation- 
dissipation theorem, D and 7 are related through the 
Einstein relation D = k B T/^. We assume that F(t) 
is time-dependent, but homogeneous. It describes ei- 
ther an externally applied force, or some intrinsic non- 
equilibrium fluctuations of the system (see below). Due 
to the single file nature of the glycerol transport through 
GlpF, a force applied at either end of the channel will 
be transmitted along the file without significant loss in 
intensity (incompressibility of the single file), which justi- 
fies our assumption that F(t) is homogeneous along the 
channel. For a periodic f(x), the above LE describes 
what is often referred to as a fluctuating force (tilting) 
ratchet [l7| . 

At this point we introduce dimensionless units that 
will be employed throughout this paper, unless otherwise 
stated. All other units can be expressed in terms of the 
following three: length of GlpF C = L « 4.8 nm, diffu- 
sion time T — tb — L 2 / D « 10 _T s, and thermal energy 
£ = k B T « 4.28 x 10" 21 J; here k B is the Boltzmann 
constant, T = 310 K is the physiological temperature, 
and D « 2.2 x 10 -10 m 2 /s is the effective diffusion coef- 



ficient of glycerol inside GlpF 16]. Thus, the force unit 
is T = jD/L = k B T/L « 0.9 pN. In the new units, 
the Fokker-Planck equation (FPE) corresponding to the 
above LE reads 

d t p(x,t) = -d x J(x,t), (la) 
J(x,t) = -d x p{x,t) + [f{x)+F{t)]p{x,t) (lb) 

where p(x, i)dx is the (unnormalized) probability of find- 
ing a glycerol molecule in (x,x + dx) (see FigQJ, and 
J(x, t) is the local, instantaneous flux of glycerol through 
the channel. The probability density p(x) is related to 
the local concentration C(x) by p(x) = S(x)C(x), with 
S(x) the area of the channel cross section. From the crys- 
tal structure one finds that the opening area at both ends 
of GlpF is So = 5(0) * 5(1) * 100 A 2 0. 
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FIG. 2: (a) Ratio of outward and inward glycerol fluxes in 
GlpF as function of a constant force Fq. (b,c) Transport in- 
duced by a square- wave force in GlpF. Shown are the inward 
[pi = 0], outward [po = 0], and equal concentration [pi = po] 
fluxes vs force amplitude Fq. (d) Ratio of the inner to outer 
glycerol concentrations vs Fq for vanishing flux. 

GlpF can be regarded as a nanopore, which connects 
two reservoirs of glycerol molecules located at x = 
(periplasm) and x = 1 (cytoplasm), respectively (Fig[Q. 
The glycerol uptake is driven by a concentration gradi- 
ent: glycerol concentration (and therefore po = p(0)) is 
finite in the periplasm and vanishingly small in the cy- 
toplasm (pi = p(l) w 0). Indeed, once glycerol enters 
the cytoplasm it gets phosphorylated by glycerol kinase 
(GK) and, as a charged particle, the product glycerol 
phosphate cannot leave the cell |2l[. However, exces- 
sive accumulation of glycerol phosphate may result in 
cell poisoning and death. To avoid this, the enzyme GK 
is genetically turned off, preventing further glycerol phos- 
phorylation [2l|. Since the glycerol concentration gradi- 
ent across the cell membrane persists, one expects that 
glycerol uptake should continue, in spite of its poten- 
tial damaging effect on the cell, until this gradient van- 
ishes (i.e., p\ = po). Below we demonstrate that channel 
asymmetry, combined with non-equilibrium fluctuations, 
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can stop glycerol uptake against the persistent concentra- flux through GlpF in four distinct cases corresponding 
tion gradient, keeping the cytoplasmic glycerol at a level to suitable choices of F(t). 
Pi < Po- To this end, we calculate the steady glycerol 




FIG. 3: (color online). Contour density plots of the numerically evaluated (a) relative inward J(Fo,To|po, 0)/Jo and (b) outward 
— J(Fo,To\0,po)/ Jo glycerol fluxes in GlpF, and (c) concentration ratio pi/po = Ao/Ai at which current reversal (J = 0) takes 
place in the channel, as a function of the amplitude Fo and mean switching time To of the RTF. Jo = J(0, Tb |po,0) is the 
inward flux in the absence of the RTF. The scales on the left and right refer to (a)-(b) and (c), respectively. 



1) Transport driven by concentration gradient [F(t) — 
0] - In the steady state, the flux is constant throughout 
the channel [c.f. Eq^H an d, i n the absence of the external 
force F(i), is given by 



J(po,Pi) = A oPo ~ Axpx, 



An = e 



Ui 



e u ^dx 







(2a) 

z = 0,l (2b) 



Since the PMF U (x) vanishes on both sides of the chan- 
nel, i.e., U = U(0) = 17(1) = Ui = (see Fig.p, one 
has Ao = A\ , and the flux is proportional to the glycerol 
concentration difference, i.e., J(po,Pi) = Aq (po — Pi) — 
A§S§ (Co — C\). Hence, the flux vanishes in the absence 
of a concentration gradient, and the flux is insensitive 
to the asymmetry of the channel. In order to produce 
directed transport, one needs to drive the system out 
of equilibrium, e.g., by applying an external force or by 
subjecting the system to non-equilibrium fluctuations. 

2) Transport driven by potential gradient [F(t) = Fo = 
const.] - The dependence of the flux on the asymmetry 
of the channel is manifest when Uq ^ E/i, i.e., in the 
presence of a potential gradient. The constant external 
force leads to an effective, tilted PMF, U e ff(x) = U(x) — 
Fox, and according to Eqs.Othe stationary flux reads 



A (F ) = 



f 

Jo 



A 1 (F ) = A (F )e- F ° 



J(F \p ,p 1 )=Ao(F )p -A 1 (F )p 1 



(3a) 



(3b) 

For a symmetric PMF, i.e., U(l — x) = U(x), follows 
Ai(-Fq) = Ao(Fo), implying that the (magnitude of the) 
flux is also symmetric | J (F \po,Pi) \ = \ J (-F \puPo) \- 
In general, for a generic asymmetric PMF the inward and 
outward fluxes will be different. Under normal conditions 
when pi = 0, the inward flux is J > = ^0(^0)^0- The 
same flux through an inverted channel, i.e., (po : pi) — > 
(0,p ) and F -> -F , would be J< = A 1 (-F )p . It 
can be readily checked that J^/Jo > 1, i.e., the con- 
stant driving force enhances the flux when applied along 
the concentration gradient. Furthermore, according to 
Fig. under identical conditions the flux through the 
inverted channel is always larger than the flux through 
the normally oriented GlpF, the ratio of the two increas- 
ing monotonically with Fq. 

3) Transport driven by an external periodic driv- 
ing force - Next, we consider an external force F(t) 
that switches periodically between ±Fq, Fo = const 
(square-wave force). Although the time average of 
F(t) is zero, this force induces a finite flux through 
GlpF even in the absence of a concentration gradi- 
ent. Indeed, assuming that one can neglect the tran- 
sient in the instantaneous flux after switching F(t), the 
mean flux J through the channel can be expressed as 
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J{Fo\Po,Pi) = [J(F \po,pi) + J(-F \po,pi)]/2, where 
J_(Fo\Po,Pi) is given by Eq. 03 Then, for p 1 = p 0: 
J I J = J(Fo\po,po)/ J(0\po, 0) is negative, and decreases 
monotonically with Fq, as shown in Fig. 03. In this 
case too, J < implies that for GlpF the outward 
flux J < = — J(— Fo|0,po) is bigger than the inward flux 
J > = J(Fo|po,0) as shown in Fig. [SJd, and both fluxes 
are bigger than Jo, the flux in the absence of the ex- 
ternal force. Furthermore, the concentration ratio pi/po 
at which the flux through the channel vanishes (current 
reversal), i.e., J(Fo\po,pi) = 0, is plotted as a function of 
Fo in Fig. [2]i. The values pi/po < 1 found are expected 
and consistent with the fact that for the same force level 
and concentration gradient, the outward flux is larger 
than the inward flux. 

4) Transport driven by non- equilibrium fluctuations - 
Finally, we consider the effect of non-equilibrium fluc- 
tuations of the cell membrane on the glycerol trans- 
port through GlpF. We model such fluctuations by a 
random telegraph force (RTF), i.e., a homogeneous di- 
chotomous force F(i), which switches between two states 
±Fo with switching times that obey a Poisson distribu- 
tion. For the RTF holds (F(t)) = 0, and (F(t)F(0)) = 
Foe -2 ^ T °, where To is the mean switching time. The 
stationary FPE in this case consists of two coupled equa- 
tions, -&lp± (x) + d x [f(x)p± (x)} ± F p± (x) - p± (a) /T + 
p T (x)/To = 0, where p±(x) is the conditional probabil- 
ity density that the RTF is in the ±To state. The cor- 
responding flux is J = —p f (x) + f(x)p(x) + FoAp(x) — 
const, where p = p+ + p- is the total probability den- 
sity, and Ap = p+ — p_. For the feature-rich GlpF 
PMF (Fig. I!]) the flux J(F 0l T \p 0lPl ) = Ao(F ,T )p - 
*4i(To,To)pi needs to be computed numerically, e.g., 
by employing a matrix continued-fraction method. The 
computed flux is shown in Fig. [HI as two dimensional 
density plots for F G [0,8] and T G [10~ 2 , 10 5 ] (log- 
arithmic scale). In SI units, these values correspond to 
T G [0,7.2] pN and T G [10~ 9 , 10~ 2 ] s, respectively. 
The conclusions drawn from Fig. |3| are consistent with 
our previous findings for externally applied deterministic 
forces. First, the RTF-induced asymmetry between the 
inward and outward fluxes is manifest (Figs. EK-b), with 
bigger outward flux for the same To , To and concentration 
gradient. Second, the RTF-induced flux enhancement is 
more pronounced for slower fluctuations and for larger 
Tq. Third, the flux in the absence of a concentration 
gradient across the membrane is always outward. Thus, 
the ratio of the inner to outer glycerol concentrations at 
which glycerol uptake ceases is always less than one, as 
shown in Fig.O^. 

Our calculations have demonstrated that non- 
equilibrium force fluctuations acting on glycerol in GlpF 
can have an important effect on the glycerol uptake by 
a cell. On the one hand, slow, large amplitude fluctua- 
tions enhance the concentration gradient-driven glycerol 
uptake (Fig. 0t), which may be vital for the cell under 



poor nutrient conditions. On the other hand, when glyc- 
erol is abundant, fluctuations provide an effective pro- 
tection mechanism to the cell by stopping glycerol up- 
take well before the cytoplasmic concentration reaches 
the periplasmic one (Fig. Oh). The mechanism underly- 
ing this behavior is related to the ratchet effect and de- 
pends sensitively on the asymmetric shape of the PMF 
characterizing glycerol conduction in GlpF. The effects 
described should be testable experimentally. In fact, re- 
cent experiments have demonstrated that cell membranes 
are subject to slow, kHz frequency (To ~ 10 -2 ) non- 
equilibrium fluctuations that may be related to pumping 
mechanisms by which cells supplement the passive diffu- 
sion of nutrients |4| . 
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